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Above-threshold ionization with highly charged ions in super-strong laser fields: 
III. Spin effects and its dependence on laser polarization 
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Spin effects in the tunneling regime of strong field ionization of hydrogenlike highly charged ions in linearly 
as well as circularly polarized laser fields are investigated. The impact of the polarization of a laser field on the 
spin effects are analyzed. Spin-resolved differential ionization rates are calculated employing the relativistic 
Coulomb-corrected strong-field approximation (SFA) developed in the previous paper of the series. Analytical 
expressions for spin asymmetries and spin flip probability, depending on the laser’s polarization, are obtained for 
the photoelectron momentum corresponding to the maximum of tunneling probability. A simpleman model is 
developed for the description of spin dynamics in tunnel-ionization, which provides an intuitive explanation for 
the spin effects. The spin flip is shown to be experimentally observable by using moderate highly charged ions 
with a charge of the order of 20 and a laser field with an intensity of / ~ 10^^ W/cm-. 


I. INTRODUCTION 

Due to advances of the laser technology, the relativistic 
regime of laser-atom interactions is now within the experimen¬ 
tal reach IfT HTOll . Recently different tools for advanced studies 
of particles in extreme laser fields have been also developed, 
which make it possible to select and prepare well-defined ion 
ensembles and to optimize the laser-particle interaction IfTTl . 
This will allow the investigation of the relativistic strong field 
ionization dynamics, which is accessible only in a combination 
of strong lasers and highly charged ions. 

The strong field approximation (SFA) IfT^ [T3ll and the 
imaginary-time method (ITM) ifTrUTSll are well developed 
theoretical tools for the analytical investigation of strong field 
ionization in the relativistic regime. In particular, the differ¬ 
ential and total ionization rates were calculated using these 
methods. The photoelectron momenrnm distribution in the 
relativistic regime of above-threshold ionization was well ex¬ 
plained by those calculations lfT9lE5l . However, the details of 
the electron spin dynamics in the relativistic tunnel-ionization 
regime still needs further elucidation. 

Spin effects in different processes in laser fields have been 
investigated since the invention of the laser |2^ . In particular, 
it was shown that the relativistic dynamics of free electrons in 
a strong laser field is disturbed by spin induced forces lIZTlI^ . 
and that the electron radiation can be modified due to the spin 
induced dynamics Il29l . Spin effects in laser assisted Mott scat¬ 
tering and laser assisted Moller scattering were investigated 
in BOIImI and 0^ . respectively. Polarization effects in the 
multiphoton Compton scattering were investigated in EMD- 
Transfer of polarization from the laser beam to positrons via 
Compton scattering and pair production is shown in ll4^ . Re¬ 
cently, a spin flip effect was shown in the Kapitza-Dirac effect 
Il43l - |45l . Relativistic spin operators in various electromag¬ 
netic environments have been discussed in Il46ll471 . Collapse 
and revival of the spin precession in a laser field has been re¬ 
vealed in ll48ll . Furthermore, interesting polarization effects 
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were explored in electron-positron pair production processes 
in ultra-strong laser fields, in particular, during electron and 
positron pair creation in combined Coulomb and strong laser 
fields ll49l . in the multiphoton Bethe-Heitler process ll50l and 
in two counterpropagating laser pulses ED. Spin correlations 
in electron-positron pair creation by a laser pulse and a proton 
beam was examined in E2l . 

During the relativistic laser-atom interaction spin effects 
were shown to appear in the laser-driven bound electron dy¬ 
namics ll53l and, in particular, in the radiation of high-order 
harmonics ESEll- The spin asymmetry in the strong-field- 
ionization process of an atom with a circularly polarized laser 
field was discussed in ESl . neglecting the spin dynamics in 
the bound state. The latter, however, can have a significant 
impact on spin effect as it was shown in E9l . The spin dynam¬ 
ics in nonsequential double ionization of helium was consid¬ 
ered in IIMIIMII . Furthermore, the photoelectron spin polariza¬ 
tion can also arise because of the electron-ion entanglement 
MM- The spin-asymmetry in the relativistic regime of 
tunnel-ionization from p-states was explored in IMl . 

In this paper we investigate the dependence of the spin ef¬ 
fects during tunnel-ionization on the polarization of the laser 
field. The spin-resolved differential ionization rates of a highly 
charged hydrogenlike ion from the ground state in a strong 
laser field of linear and circular polarizations are calculated 
using the relativistic Coulomb-corrected SFA (cc-SFA) devel¬ 
oped in paper II (the second paper of this series ll^l66l ). Spin 
asymmetries and spin-flip effect during direct ionization of an 
hydrogenlike system are investigated for the peak of the final 
momentum distribution. Similarities and differences of the 
spin asymmetries and spin-flip effect in the cases of linear and 
circular polarizations of a laser field are analyzed. In addition 
to the standard relativistic Coulomb corrected SFA (s-cc-SFA), 
we also apply a dressed cc-SFA (d-cc-SFA), which is based on 
the use of a non-standard partition of the Hamiltonian within 
the SFA formalism lICTl l68l . The physical relevance of the 
different versions of the SFA formalism is discussed. While in 
s-cc-SFA the influence of the laser field on the electron spin 
evolution in the bound state is not taken into account, it is 
fully accounted for in d-cc-SFA, which is shown to have a 
decisive impact on the spin effects. Finally, we provide a sim¬ 
pleman model for intuitive understanding of the spin effects. It 
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incorporates the propagators of the spin states for the bound 
and continuum motion, and quasiclassical description of the 
tunneling process. 

Spin effects in the tunneling regime of ionization emerge 
through three steps ll59l : spin precession in the bound state, 
spin rotation during tunneling, and spin precession during 
the electron motion in the continuum. In s-cc-SFA the spin 
dynamics in the bound state is completely neglected. Therefore, 
in this case the spin effects are determined by the electron 
dynamics during tunneling and the motion in the continuum. 
Because of the evident asymmetry in the spin evolution in 
this picture, relatively large spin effects arise. However, it is 
known that the laser field can induce a large spin precession in 
the bound state Il54ll55]l . which is accounted for in d-cc-SFA 
t5^l6^ . It reduces the asymmetry in the spin dynamics during 
ionization and, consequently, lead to a reduction of spin effects. 
Analysis based on our simpleman model shows that the spin 
asymmetries are a consequence of the tunneling step. 

The plan of the paper is the following. The spin-resolved 
differential ionization rates within the two versions of cc-SFA 
(standard and dressed) are calculated in Sec.[^ The bound state 
spin dynamics, which is essential for d-cc-SFA is investigated 
in Sec. nni The final momentum distribution of the tunnel- 
ionized electron is presented in Sec. IV which later is used in 
the calculation of the spin effects for the maximal mnneling 
probability. Analytical formulas for the spin asymmetries and 
the spin flip in linearly and circularly polarized laser fields are 


calculated in Sec.|V 


A simpleman model for the spin dynamics 


is presented in Sec VI The possibilities for an experimental 
observation of spin effects are discussed in Sec. VII Our con¬ 
clusion is given in Sec. VIII Atomic units (a.u.) and the metric 
convention g - (-I-, -, -, -) are used throughout the paper. 


II. THE COULOMB CORRECTED STRONG FIELD 
APPROXIMATION 


The Hamiltonian which governs the dynamics of the laser- 
induced tunnel-ionization from the ground state of a hydrogen¬ 
like ion, is given by 

H^ca-(p + Ai)-Al+l3c^ + V{r), (1) 

where the laser field is described by gauge potentials given in 
Goppert-Mayer gauge as = (A°, cA/,) = -x ■ 
with the wave vector = w/c(l, k), the laser frequency w and 
the phase 77 = kx/w, and the Coulomb potential V{r) - -xjr, 
K being the charge of the hydrogenlike ion, and /?, a are the 
Dirac matrices. 

The transition amplitude between the initial state | with 
the magnetic spin quantum number s and the final state | ) 

with the number s' can be written as 

X OO 

( 2 ) 

00 

with the interaction Hamiltonian Hi„,. The S-matrix treatment 
is exact as far as it incorporates the exact final state | 'Ej'), which 


is the exact solution of the Schrddinger equation 

(3) 

The initial state in its turn fulfills the following equation 

idr\^^) ^ (H - . (4) 

In the SFA, the exact final state | 'Pj') in the transition ampli¬ 
tude 0 is approximated by the Volkov state ll69ll . whose wave 
function in the Goppert-Mayer gauge reads 


Pv — \i~ sxp {iS + iA ■ x) 
'' e 


1 + -— (1 -H k ■ or) A ■ a 
2cA 


with the Volkov action 


s in) = -px - T 

A 


uy 


A{ri') ■ p + 


Aiif 


Vs' , 

(5) 

( 6 ) 


and the free particle spinor 


V.v' = 


e + 
2 c 2 




-^p ■ o-x. 


(7) 


Here a = sjc'^ + c^p^ is the electron energy, A = sjc^—p-kjc is 
the integral of motion for the electron in a plane wave field, A s 
- J^^E(jj')dri', and x+ - (1 0 )^ and;^^- = (0 1 )^ are the 

two components spinors. Note that the Volkov wave function 
in the Goppert-Mayer gauge is obtained by first solving the 
Dirac equation in the velocity gauge A^{ri) = (0, cA{ri)) and 
then applying a gauge transformation with the gauge function 
A{t]) ■ X. 

In this approximation, the transition amplitude neglects the 
effect of the Coulomb potential on the electron dynamics in 
the continuum as well as the influence of the laser field on the 
bound state dynamics. To account for the Coulomb potential 
on the electron dynamics in the continuum, the relativistic cc- 
SFA was developed in paper II ll 66 l . Rather than the Volkov 
solution for the continuum electron Py , cc-SFA employs the 
wave function of the electron in the laser and Coulomb fields 
in the eikonal approximation 11701 - 17411 . which is given by [see 
Eq. (11.29)] 

P^ ^P^exp[iS,(x,,])], ( 8 ) 

where 

r°° B(n') 

Sc(x,J]) = j dj]'-^V{x(T]')), (9) 

with the relativistic trajectory of the electron in the laser field 
x(ri') - X + fj di]” p{tj") I A, and the energy-momentum of the 
electron in the laser field 


Pid) ^P + Mif) + k 
e(rj) - a + 


[p + A{^i)|2\■A{l^) 


C/1 

[p + A(?7)/2] ■ A(77) 


( 10 ) 

( 11 ) 


A 
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where we have used Eqs. ( |205| l and ( |206| l of Appendix |IX A| and 
define the final values of the physical variables as p = pipf), 
e = sirif), and Aipf) = 0 . 

Thus, the transition amplitude for the strong field ionization 
in cc-SFA reads 

X OO 

( 12 ) 

oo 

Accordingly, the spin resolved differential ionization rate for a 
certain spin transition can be defined as 

(13) 

d^p n 

where the rate is averaged over a laser half-cycle. 

In contrast to the exact S-matrix treatment, the results of the 
SFA calculation depend on the partition of the full Hamiltonian 
lEl, i.e., on the identification of the interaction Hamiltonian 
///„,. In the next sections we specify two different choices of in¬ 
teraction Hamiltonians, which yield different SFA approaches. 
We will calculate the spin resolved differential ionization rates 
in these two approaches and will later discuss their physical 
relevance. 

We calculate the ionization rates in linearly as well as circu¬ 
larly polarized laser fields. In the velocity gauge, the vector 
potential of the laser field can be written as 

Eq 

A(t]) = — [sin(w 77 )jc - ^ cos(w 7 ;)y] , (14) 

U) 

which yields the corresponding electric and magnetic fields 

E{ri) = -Eq [cos(m 77 )^ -H ^ sin(m 77 )_p] , (15) 

B{ri) = Eq sin(ciJT/)x - cos(m 77 )y] , (16) 

respectively, with the laser field amplitude Eq and the polar¬ 
ization parameter such that f = 0 corresponds to linear and 
^ = 1 is for circular polarization of the laser field. Further, 
we specify the propagation direction ask - z, which implies 
7] - t- zlc. 


A. Standard cc-SFA 

In s-cc-SFA, the total Hamiltonian is partitioned as follows 

H^H^q + HI,, (17) 

Hq - ca ■ p + fic^ + V(r ), (18) 

Hl,^x-E[\-a-k) . (19) 

In this partition, the initial state, fulfilling Fq. 0 . is the ground 
state of a hydrogenlike ion whose position representation 
is 

AH I 2-I„lc^ 

tlf-Qix) = -p -v — ——j-(2Kr) Ed- exp {-xr - isot) u ,, 

\ r{3 -2IplA) 

( 20 ) 


with the ground state spinor 


A* 


l—x ■ (TXs 
'' CK 


( 21 ) 


the ground state energy bq - — Ip, and the ionization energy 

Ip — ll75l . Then, after plugging the eikonal- 

Volkov ([^ as well as the ground state ( [20| wave functions into 
the transition amplitude, and changing the variables from (f, x) 
to ( 77 , Jc), we obtain 


^s-^s 

' = A 


d^xe 


I dT]e‘^’-’^'> (22) 

\J ~00 

yt (i _ q. . k)iH , 


with the following contracted action S {rf) and the relativistic 
kinetic momentum ^( 77 ); 


5 ( 77 ) = ^ J" drj' ^A(77') ■ p + + 4(s - eo)j , (23) 

qirf) ^ p + A{t]) - - — —k. (24) 

c 

Here, the prefactor is 


N = -7 A — 


C 2 


Iplc^ 


s \ r(3 - 2Iplc^) 


( 2 k) 


rhid 


(25) 


For the later convenience, we can also write down the con¬ 
tracted action in the form of 

5 ( 77 ) = ^ dr]' [q^ii]') + K^\ . (26) 

The 77 integral in Fq. (22 1 is performed via the saddle-point 
approximation (SPA). The saddle point equation, S ( 77 ^) = 0, 
yields 


q ins) + 7 f^ = 0 . 


(27) 


As in the SPA-contour only the integration region near the 
saddle point makes the main contribution to the integral, the 
contracted action can be expanded around the saddle point 77 ^ 
as 


5(77) = 5(77,) + 5(77,)(77 - 77,)2/2. (28) 

Then, we replace 77 with 775 in the rest of the phase dependent 
functions, except in the term qirf), which after the coordinate 
integration yields to a singular function at the saddle point, 
as we will see below. Afterwards, the transition amplitude 
becomes 


X CXD 

dq exp [iS ( 77 ,) -H i§ (j],)(j] - 77 ,)^/ 2 ] 

OO 

-im-x-xr+iSAx,n,)p-iAc\ . Eiq,)v\(l - a 


^s—>s 

- = A 


d^xe 


(29) 
k)us , 
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where the Coulomb correction factor arises [see Eq. (11.57)] 
Qr = exp[i5c(-»:,77,)] 




\ c 


X ■ E{tJs) 
4/„ 




Furthermore, using the fact that x ■ E/r ~ 1, see 
transition amplitude reads, 


(30) 
, the 


M,_, = <31) 


giS(ris)+iS(ris)(ri-risf II 


where 


|£(^.)l 


iisiri) = 


-/o(?7)vl, (l -a-k)us{Tf), 


Xs 

.Ip J\{rf) ■ (T 

I - TT^Xs 

. CK Jo(T]) 


with 


Jo(Jl) - J" d^x exp [-iq(T]) ■ x - Kr\ , 
Jiirj)- (T - J"(Px exp [-iq{r]) ■ x - Kr]x ■ (r. 


(32) 


(33) 

(34) 


With the help of the plane wave expansion, these space integrals 
can be calculated exactly as 


Join) 


Sttk 


(35) 


[qiqP + K^f 
i%n 

J\in)-(x^-^—^ - —2qiq)o-. (36) 


[qiqp + 

Then, the transition amplitude reads 


where we have used qirjs) - ix. Finally, the transition ampli¬ 
tude can be written as 


M, 


exp(/ 5 (? 7 .)) ^ . 

[q(7j,) ■ Vd ^ 


-a ■ Usiqs), 
(40) 


with N = iN{2mpl^Ke^’pl‘J (l - ^) 

In order to evaluate the transition amplitude for any spin 
quantization axis, we can use the rotation operator G after 
fixing the representation. Namely, we can choose the repre¬ 
sentation of the gamma matrices in the z-basis and specify 
the two-component spinor in Eq.(|^ as well as in Eq.(|21|l 
along the z-direction. Then, the spin states along an arbitrary 
quantization axis can be written as 

s' 

where is the state whose spin quantization direction is 
the z-axis and Gss’iS, (f>) is the Wigner D-matrix, which can be 
defined as 




' cos(|) e ''*sin(|)' 


(42) 


,e''^sin(f) -cos(f), 

with the spherical coordinates 0 and 0 , see Fig. [T] In other 
words, the rotated states expressed in terms of the z basis read 


r)^cosQ 

|'F;) + e'^sin(^ 


(43a) 

W-) = e^''^sin| 



(43b) 


As a result, the transition amplitude valid for any spin quanti¬ 
zation axis can be written as 


X “ giS{ri,)+iS(-n,)(ri-r],fl2\^£f^P^s^]Iplc 

... 

cx 


dq ■ 


[q{qp + Af2] 


v[, (l -a-k)u,{q). 

(37) 


As a final step, we can evaluate the 77 -integral. The pre¬ 
exponential integrand has a singularity at the saddle point 
because of the saddle point condition Eq. (27 1 . Therefore, we 
first expand the singular factor around the saddle point as 


(44) 

ij 


with being the transition amplitude when the spin quan¬ 
tization axis is chosen along the z-axis. Furthermore, the spin 
resolved differential ionization rate for an arbitrary spin quanti¬ 
zation axis can be found via 


dW.,^s' 

d^p 


2 

yij,kj 


(45) 


[qinP^ + x'^f 4 [^( 77 ,) ■ 4 ^( 77 ,)]^ (77 - 77 j )2 

and include it in the integration following the modified SPA 
1761. The bispinor Usiq) has no singularity at the saddle point 
and reads 


Usiq.,) 


Xs 


■ h 

i — 
CK 


q(0s) ■ o-xs^ 


(39) 


B. Dressed cc-SFA 

Although, s-cc-SFA improves the results in comparison to 
the usual SFA via the Coulomb correction to the continuum 
electron wave function, it neglects the influence of the laser 
field on the bound state. Because of that, the electron spin 
dynamics in the bound state is completely neglected and the 
electron spin in the tunneling bound state is the same as in 
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FIG. 1. Spin states along an arbitrary quantization axis can be ob¬ 
tained by acting the rotation operator on the initial spin states whose 
quantization axis is along z. The rotation is defined by angles (6, ^). 
The configuration of the laser fields at the instant of ionization is 
specified. 



the initial state before the interaction with the laser field. In 
this case the spin effects are determined solely by the electron 
dynamics during the tunneling and during the motion in the 
continuum. Due to the evident asymmetry in the spin evolution 
in this picture (frozen spin in the bound state, rotating spin in 
the tunneling, and oscillating spin in the continuum) relatively 
large spin effects arise. On the other hand, it is well-known 
that the laser field can induce a significant spin dynamics in 
the bound state l54ll . Moreover, the Zeeman-splitting of bound 
state levels can have an impact on the tunneling probabilities, 
in this way modifying the spin effects Czl. Therefore, it is 
important to take into account the laser field influence on the 
spin evolution in the bound state when calculating spin effects 
in the ionization process. 

With this motivation, we employ d-cc-SFA, see paper II, 
which is based on a specific partition of the Hamiltonian, in 
which the bound state dynamics in the laser field is accounted 
for. In d-cc-SFA the total Hamiltonian is split up as follows 

(46) 

Hq - ca ■ p+ + V(r) - {x ■ E){a ■ k) , (47) 

Hl^x-E. (48) 

In this case the SFA transition amplitude reads 

dtJd^x'Fl.'\x,t)(x-E)0^g(x,O, (49) 

where the so-called dressed ground state satisfies the following 
equation 

= (50) 

This equation includes the spin precession in the bound state 
induced by the laser field. 


The Schrodinger equation for the dressed bound state (50 1 
cannot be solved analytically in an exact way and, therefore, 
several approximations are applied. Since the typical dimen¬ 
sion of the atomic bound state is much smaller than the wave¬ 
length of the laser, we apply the dipole approximation for 
treating the dressed bound state dynamics, i.e., rj ^ t. Further¬ 
more, we are concerned mostly by the electron spin precession 
in the bound state, consequently, we consider only transitions 
in the subspace of spin states and describe the dressed ground 
state with the following ansatz 


10^(0) = 2 (Olios'(0), 


(51) 


where |i/^g (f)) is the ground state wave function whose position 
representation is given by Eq. ( |20| l. According to the wave 
equation (501, the coefficients (f) satisfy the following dif¬ 
ferential equation 


/C" 


"(t) = 2 C^^'(f)<<(0|//i(OI(Ao (0>, (52) 


where Hpt) = -E{t) ■ xa ■ k. When the spin quantization 
direction is the z-axis, the value of the matrix elements can be 
written as 


<-Ai(0|//i(f)llAj(0) = <-Ao(0|//i(OI(^o(0) = 0, (53a) 

= . lEyio-EAt) ^^ 

^ (53b) 


where b s 1 -2Ipf{3c^). Thus, we obtain the following coupled 
differential equations 

( 7 +±(f) = , (54a) 

C-^(t) = C-^(t)F^(t ), (54b) 


with 


^ iEyjt) ± EAt) 
^ ’ 2c16 


(55) 


It is useful to separate the bound state propagation in two 
stages; first, from the switching on the laser field up to the ion¬ 
ization moment, and second, the under-the-barrier dynamics. 
The description of the first stage requires the exact solution 
of Eq. ( |54| l, which for linear and circular polarizations of the 
laser field are discussed in Sec|III| Eor the second stage of 
the under-the-barrier dynamics during the imaginary time, a 
simple approximate solution of the bound state wave function 
can be found exploiting the shortness of this time propaga¬ 
tion. In fact, the time integration in Eq. ( [49] l will be performed 
by SPA and the matrix element ( |49| ) will be evaluated at the 
saddle point tJs - rjr + irji (see Sec |IV[ |, where rjr is the in¬ 
stant of ionization at which the spin state leaves the bound 
state, and 77 , is the duration of the imaginary time (the Keldysh 
time), during which the state evolves under-the-barrier, yet, it 
is small with respect to the laser period in the tunneling regime 
\u)rii\ ~ j <s: 1 [the tunneling regime of ionization is defined by 
the Keldysh parameter y = yfUpColEo « 1]. 
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In order to investigate solely the under-the-barrier propaga¬ 
tion, we first convert these differential equations to Volterra 
integral equations of the second kind, which for the time prop¬ 
agation after f = can be written as 

C+*(f) = -H C,f J ds F^{s) (56a) 

+ J dsC^^(s)F^(s)J^ dTF^(T), 
c~\t) = J dsF^{s) (56b) 

+ J dsC-^(s)F^(s)J^ drF^Cr), 

where Cg^ s (tr). For calculation of the maximal tunnel¬ 
ing probability, we choose = 0 , which corresponds to the 
electric field maximum, see Eq. Then, introducing the 
dimensionless parameter ip = tot, we arrive at the following 
integral equation 


{tp) = r iv) + 




duK’^^ {(p, m)C” (m) , (57) 


with the relativistic invariant field parameter ^ s EoKcco). 
Here, the integral Kernels are given by 

u) - K^~*{ip, u) - [cos(m) + i( sin(M)] (58) 

X {i( [cos(m) - cos(i/j)] + sin(M) - sin(i/j)) , 

and the functions /” (f) are 

r^(^) - + Cr f [*^(cos(^) - 1) + sin(^)] , (59a) 

/-*(^) - Co* + Co^^ [i^(cos(^) - 1) + sin(^)] . (59b) 

The formal solution to the integral equation ( [57] ) can be 
given by the iterative method iTTSlI . Although in weak laser 
fields, ^ « 1, a few iterations would give an accurate result, 
for the strong fields, ^ » 1, generally, one has to deal with 
the infinite sum of iterations. Nevertheless, we consider a very 
short time propagation <p — cotji « 1. Moreover, we observe 
from Eq. that the effective time parameter scales as 

^6(p, which can be estimated as ~ ^^pjc 1 and it is 
small. Therefore, the solution of the integral equation ( |57| l 
for the considered small time propagation can be represented 
by the leading iteration term even for strong fields. The first 
iteration, which leads already to terms of the order of (/p/c^)^'^^, 
is given by the following expressions 

C**(^) = C** + i 

+ ^ [4Co*"d^ - 2iC;^d^e^ + Co*"d'^'] , 

C^\p) - Co* + ^-d^p - i [2/Co"5^^ + Co*b'f"] 

^ ^ - 2 iCo-*d 2 ^ 2 ^ + Co"d'^'] p^ . 


Furthermore, the coefficients (60 1 can be represented as 

a^'(p)^n(p,o)a/, (61) 

with the transformation matrix 


n(p,0) 


,2 \ 


1 - Id^^^p^ ^^d^p - ^^idC^p- 


24" 


-\d^p-\idc^p^ 


1 - id^ew' 




(62) 

In this way we separate the propagation of the spin states 
through the imaginary time axis from the states at the instant 
of ionization. 

Then following the same procedure as in the case of s-cc- 
SFA, the transition amplitude in d-cc-SFA for the case when 
the spin quantization direction is the z-direction can be written 
as 


Ml = A 


exp 


[iS (77s)) 


|£(77s)I 


ipM 


[q{T],) ■ Vt 


X 


1 H- —or-A(77 s)(1 -I-O' -A:) 
2cA 


2^C^0 Ul,.(qp), 
(63) 


where 


t/^(77s) 


1 - y^ficoqsf ± ^id^(f{coq,f 


Ulills) 


1 


-d^{uT]s) ± -id^^(urjs) - I I (ojfjsy 


d^f d^ 


M^T/s) . 


The transition amplitude as well as the differential ionization 
rate for any spin quantization axis can be defined via Eqs. 
and (|45]l, respectively. 


III. SPIN DYNAMICS IN THE BOUND STATE 

The spin resolved ionization amplitudes in d-cc-SFA, given 
by Eq. (63 i, depend on the coefficients Cg*, which describe 


the spin precession in the bound state when the quantization 
axis is along the laser propagation direction and are evaluated 



(60a) 

at the ionization time f^. These coefficients are found from the 
solution of the system of differential equations (54 1 with the 
following boundary conditions 



limC**(f)=l, (64a) 



limC*^(f) = 0, (64b) 


(60b) 



describing the initial spin state of the atom when the laser field 
is switched off adiabatically. 
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A. Circular polarization 


In the case of circular polarization, the solution of Eq. ( |54| ), 
which satisfy the boundary conditions (64 1 , is 


(0 = 


+ (i + vi+ 

i5^ 

^(52^2 + (l + Vl + 


(65a) 


± I 1 + 




(65b) 


The maximum of the ionization probability is at the instant of 
ionization = 0, see Eq. (0, at which the coefficients yield 


c:utr - 0) - C-- (t, = 0) = 


c^-^r = 0) = c-.;,(f, = 0) = 


1 + 


+ (1 + Vl + 

(66a) 

__ 

yjs^f + (1 + Vl + 

(66b) 


In the weak field limit, ^ « 1, one has the following asymp¬ 
totics 


c++(f, = 0) = C--(f, = 0) = 1 - yv , f « 1, (67a) 

C+-(f, = 0) = C-^(tr = 0) = ^ « 1, 

I io 

(67b) 


In the weak field regime the averaging result coincides with 
the instantaneous one with an accuracy of the order of TfTo 
(TITq <Si 1). However, in the strong field regime the period of 
the spin oscillations are ^ times smaller than the laser period. 
In this case we average the ionization probability over the time 
period of the spin oscillation T - TqI^. 

Erom Eq. ( |63| l one can see that the spin dependent part of 
the transition amplitude in the case of d-cc-SEA is determined 
by the following factor 

Y, C^^"(tr)U]M ), (69) 

S" 


with 


y]Hri's) = 


,2t 


1 H- —a-A(r]'^)(l + ak) 
2cA 


(70) 


Here, we did the replacement Cq " C^^"(tr) as well as 

tr + T]s in order to investigate effect of an arbitrary instant of 
ionization tr. Then, the averaged differential ionization rate 
contains the factor 

= \SUsW\\\Cltrcitrt) + 

+ {i,){C:l,{tr)C-:{tr)) 

+ , (71) 

with = 12y^(77j)t/j(7/'). Here, the slow oscillating 

functions *S^_^^,(7s) with respect to the averaging time T are 
taken out from the averaging. Erom Eq. ( |65| l one can see that 
the mean values are time independent, and the fre¬ 

quency of the oscillations of C^J'^^(tr)C^J^J‘(tr) is to. Therefore, 
the averaged probability over the time T <s: Tq will coincide 
with the instantaneous value 


while in the strong field regime, ^ » 1 



c:;rcitr = 0) = C-- (C = 0) = ^ , 

V2 

f»l. 

(68a) 

c:-,(tr = 0) = c-;,(c = 0) = ^, 

V2 

f»l. 

(68b) 


We observe from the solution ( |65| l that the spin oscillates 
over time with the frequency w (l + /2. While for 

weak fields it oscillates around the value at = 0 with the 
laser’s frequency o), for strong fields the spin state does full os¬ 
cillation between the up and down states with a large frequency 
d^w/2. As the coefficients C^J^^{tr) and, accordingly, the spin 
resolved ionization probability oscillate with respect to the ion¬ 
ization time tr, it is physically more appropriate to consider the 
time averaged ionization rate. The latter, in fact, corresponds to 
the averaging over the photoelectron momentum, because there 
is a mapping between the final momentum and the instant of 
ionization via the saddle point equation. Choosing the interval 
of the time averaging T much smaller than the laser period 
To, one can still relate the averaged ionization probability to 
the maximum of the momentum distribution of the tunneling 
electron. 


= (72) 

However, this particular choice of the spin quantization axis 
along the laser propagation direction is very special in the cir¬ 
cular polarization case. When the quantization axis is arbitrary, 
the probability is a linear combination of (C)C'*;rc^*(C) 

with arbitrary values of i, because 

i,j,kj 

X {Sl^j{ri:)CrUtr) + S^_^j(rj:)C-r,{tr)) 
which, for instance, contains a term 

<si_//7.)^r-./(7.)<Cri;c(f.)C':;(c)) ■ ( 73 ) 

The term C^~^(tr)C~^y^*(tr) oscillates with the frequency 
^6a) for strong fields, and its mean value is vanishing 
^tr)) - 0 when the averaging period fulfills the 

condition 

1 

-« r « To . (74) 
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Note that that the instantaneous value of this term at = 0 is 
1 /2 in the strong field limit, ^ » 1. Therefore, the averaged 
probability will differ from the instantaneous value in those 
cases when the direction of quantization axis is chosen other 
than in the laser propagation direction. 


B. Linear polarization 


In the case of linear polarization of the laser field, the so¬ 
lution of Eq. ( [54| ), with the boundary conditions according to 
Eq. ( |64| ), is 


Cun(t) = COS 


-6^ sin(wf) 




-6^ sin(wt) 


(75a) 

(75b) 


The spin states at the instant of ionization, L = 0 are 


Cunitr = 0) - 1, (76a) 

C,^:(L = 0) = 0. (76b) 


Similar to the circular polarization case, in weak fields ^ «c 
1 , the spin states oscillate over time around the value at the 
instant of ionization L = 0 with the laser’s frequency, 

Cto(0= l-^(5^sin(cuf)/2)^ ^«1, (77a) 

Cf;:(t)^+y^sm(ajt), ^«1, (77b) 


IV. FINAL MOMENTUM DISTRIBUTION OF THE 
TUNNEL-IONIZED ELECTRON 


We are concerned with the spin resolved ionization prob¬ 
abilities corresponding to the maximum of the momentum 
distribution of photoelectrons. In this section we derive the 
momentum distribution of the directly ionized electrons, taking 
into account the relativistic corrections during the under-the- 
barrier motion (~ Ipjc^) as well as the nonadiabatic corrections 
(~ y^). The momentum corresponding to the maximum of this 
distribution will be used in the next section for the evaluation 
of the spin asymmetries as well as for the spin flip. 

In SEA the dominating part of the ionization probability is 
given by the tunneling exponent 

W ~ |exp (i5 (p, tjAp)))\ , (79) 

where the exponent depends on the final momentum p as well 
as the saddle point rjs - rjr + irji. Furthermore, via the saddle 
point equation qipsf' - -k^, the momentum and the saddle 
point are connected to each other so that the exponent is a 
function of the final momentum. In other words, each saddle 
pints correspond to different momenta for the tunnel-ionized 
electron. In order to find the momentum that maximizes the ion¬ 
ization probability, we first note that in the quasistatic tunnel- 
ionization regime (y « 1 ) the tunneling probability is maximal 
when the electric field reaches the maximum at the real part of 
the saddle time, i.e., when 


\E{rir)\ = Eq . 


(80) 


Consequently, in weak fields ^ ^ 1 the time averaged probabil¬ 
ity, carried out over the time region given by Eq. ( |74| ), coincides 
with the instantaneous one. 

In strong fields ^ » 1, the spin oscillates between the up 
and down states with the frequency i5^w/2. The mean value 
of the differential ionization rate for an arbitrary spin quanti¬ 
zation axis in the linear polarization case includes the follow¬ 
ing three factors; (cos |^ 5 d^sin(cufr)j sin |^j5^sin(wL)j) = 0 , 

(cos^ |^j5^sin((yfr)j) = 1 / 2 , and (sin^ j5^sin(wL)j) = 1 / 2 , 
when the averaging is carried out over the time region ( [74) i. Ac¬ 
cordingly, when the quantization axis is along the z-direction, 
the averaged spin resolved ionization probability in the strong 
field limit is determined by the factor 

= \ f » 1, (78) 

which does not coincide with the instantaneous value of the 
transition probability. The rate for an arbitrary quantization 
axis can be calculated in a similar way. 

As a summary, the spin resolved differential ionization rate 
corresponding to the maximum of the photoelectron momen¬ 
tum distribution is well defined in the case of weak fields and 
corresponds to the instantaneous ionization rate at the peak 
of the laser field. In strong laser fields, the spin resolved ion¬ 
ization rate is highly oscillating with respect to the ionization 
time, therefore, it is physically relevant to average the rate over 
a period fulfilling the condition ([74|. 


Here corresponds to the instant of ionization in the quasi- 
classical description of tunneling such that the electron leaves 
the bound state and starts the continuum motion. Then, taking 
into account Eqs. ([T5]l and ([T6]l we find rjr = 0 , and 


E(T]r)^-EQX, (81) 

B(pr) = -Eoy , (82) 

Mpr)^ -—0 ■ (83) 


Note that for a circularly polarized field, at any time the condi¬ 
tion ( |80| ) is fulfilled. Nevertheless, in order to generalize the 
result to an arbitrary polarization, we set - 0. This condition 
further implies that the most probable tunneling is along the 
direction that points the maximal electric field. Within this 
conclusion, the saddle point equation (|27]i reads 


2^ 2/, J ,/2^2^/ £o^cosh(w77,)\^ 

K +C [I-A-Ipjc ) +lpy - - - I 


+ \Px + l 


Eq sinh(w)7,) 


U) 


= 0 . 


(84) 


For a real rji, Eq. (j 8 ^ can only be fulfilled if 


Px = 0, 


(85) 


i.e., the final momentum along the tunneling direction of the 
tunnel-ionized electron should vanish. Correspondingly, the 
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contracted action 


can be written as 


Sirjs) = 


- 1) [sin(2w77j) - cos(2w77J] 
8/iw^ 

^Eopy^to [sin(w?7^) - ojtJs cos(w77j)] 

■ 


( 86 ) 


Further, it is more convenient to maximize the tunneling proba¬ 
bility not by the final momentum p = p{rj f) but by the momen¬ 
tum pirjr) at the tunnel exit, assuming the electron propagation 
from the tunnel exit to the detector is governed by the laser 
field only. The latter is defined via the following relation. 


P = - A{T]r) - — {piT]r) ' 


Mrir) 


A(77,), (87) 


where we have used Eq. (205 i in Appendix |IX A and the fact 
that Aijjf) - 0. Taking into account Eqs. (83 i and (85 I, the 
relation ([87]l yields 


Px = PxiPr) = 0 , 


(88a) 

E En _ 

Py = Pyirir) + '- , 

OJ 


(88b) 

.X 1 ^£0 

/ , s <Eo\ 

(88c) 

Pz - PziPr) + , 

Ac CO 



We observe that the initial momenrnm along the mnneling 
direction at the tunnel exit, px{ijr), vanishes, which agrees with 
the simpleman model as rjr is the turning point, i.e., the tunnel 
exit point, where the classical particle has a vanishing velocity. 

Next we express the constant of motion A via {pyirjr), Pz(t]r)) 
as 


^ ( ^]c^ + (Pyinrf + Pzi-Tlrf) “ Fz('7r)j , 

Pzi^r) 


1 - 


(89) 


where we neglect Pyirjrf'lc^ and PziPrf 'which are of the 
orders of {EolEa){2Iplc^) [the latter follows from the fact that 
the width of the transverse momentum distribution at tunnel- 
ionization is V^/(2/p)*^"^ lITfil . Note that in the mnneling 
regime EqIE^ < 1/10, with the atomic field strength £„ = 
(2/p)^^^]. Using the A value, the contracted action ( 861 reads 


SiPs) 


Eq [Eo + C<^Py{rir)\ 
3[1-p,(77,)/c] 


7?^ 


(90) 


where we expand the action ( [86| ) up to nonvanishing order in 
rjs taking into account that \toris\ ~ y 1. At the same order, 
the saddle point is calculated as 




Py(rjr)^ + Pz(rjr)^ + 2/p[l - Pz(>?r)/c] 

' Eo[Eo + {Py(TJr)co] 


(91) 


Einally, we can identify the momentum which maximizes the 
tunneling probability via the condition 


dS 


dpjjlr) 


= 0, 


(92) 


0.1 F 


0.05 


0.00 - 


-0.05 


-0.1 b 



-0.01 -0.005 0.000 


0.005 


0.01 


FIG. 2. (Color online) Ionization probability calculated numerically 
via Eq. (|79|l versus the final transverse momentum. The final mo¬ 
mentum is defined as = Pj^(l -i- <5p^), with given by Eq. 1 95b i, 
and Eq. j95c| and being the dimensionless deviation. The applied 
parameters are k = 50, w = 1, /" = 0.5, EolEa = 1/30. Note that we 
set ai = 1 instead of ai = 0.05, which is a typical parameter for tunnel- 
ionization, in order to confirm that the derived analytical results are 
valid in a nonadiabatic regime as well. 


with Pj_(Tjr) = (PyiPr), PziPr))- The Condition yields the follow¬ 
ing two equations 

PyiPrf - 2pz(T]rf - Ip + CPziPr) (S -H Ip/c'^) = 0 , (93a) 

PziPr) 


6EoPy(rjr) + / Spyirjrf - PziPrf 


• 2 /„ 


= 0 . 
(93b) 


Up to the leading order in Ipicr, the solution of Eq. (93 i can 
be given by 


Pzi^r) = 
PyiPr) = 


jL 
3c ’ 
^tolp 

3Eo 


(94a) 

(94b) 


As a summary, the final momentum yielding the maximal 


tunneling probability, with Eq. (88 i, can be written as 


= 0, 


(95a) 



(95b) 



(95c) 

3c 2c \ 

3c^j 


the validity of which is illustrated in Fig.|^ where the final 
momentum distribution is calculated numerically via Eq. ©■ 
We emphasize that the momentum distribution (95 I is valid 


also for an arbitrarily elliptical polarization with an ellipticity 
0 < ^< 1 . 

Furthermore, by using the final momentum ( [95| ), we can 
write down the saddle point, which maximizes the mnneling 
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probability, as 


Vs = 


Eo 


5/„ 




is the total ionization rate averaged by the initial and summed 
over the final spin states. 

Secondly, one may define the spin polarization asymmetry 
parameter as 


where /p/c^-expansion is applied. Accordingly, the maximal 
probability is derived from Eq. (|79|) can be found as 


W ~ exp 


[iS (/?i))| 


exp - 


'3£o 


1 - 




(97) 

Here we should note that the initial transverse momentum 
given by Eq. ( [94] i arises during the classically forbidden under- 
the-barrier dynamics in tunnel-ionization. The momentum 
along the laser propagation direction Ipf{3c) is a relativistic 
signature of the process EU, whereas the momentum along 
the direction of the laser’s magnetic field at the instant of ion¬ 
ization ^£’oT^/(6w), which depends on the laser polarization, 
is induced due to the nonadiabaticity of the ionization and it is 
significant at large Keldysh parameters y IfSOllSTl . In addition 
to the SEA prediction, there exist also an initial nonvanishing 
momentum along the tunneling direction for the most probable 
trajectory, if one goes beyond the quasiclassical description of 
tunneling and defines a tunneling time delay Il82l . 

In the following sections, we discuss the spin dynamics in 
the tunneling regime, y I, consequently, we will omit the 
y^ terms in Eqs. (95 i - (97 1 in the corresponding calculations. 


Jlp - 


Wt 


( 101 ) 


which is a measure of the electron final polarization in the case 
of initially unpolarized states ll58l . 

In addition to the asymmetries we will provide also the spin 
flip rate 


n = ■ ( 102 ) 

Wt 

The three independent parameters and fully de¬ 

scribe the spin transitions, taking into account the normaliza¬ 
tion condition (l/2)2s,.(' Ws^s'l^^T - 1- 

To sum up, we have introduced the following two sets of 
parameters: the laser polarization parameter and the spin 
quantization angles 0, (p in order to investigate analytical results 
for the spin effects. The latter are described by three indepen¬ 
dent parameters; the tunneling and polarization asymmetry 
parameters, and the spin flip, for which we discuss the weak 
field as well as the strong field limits in the following section. 


V. SPIN ASYMMETRIES AND SPIN ELIP 


A. The s-cc-SEA prediction 


In order to investigate the spin dynamics we consider the 
following physically relevant choices of the spin quantization 
axis; along the laser propagation direction, along the direction 
of the laser electric field, and along the laser magnetic field at 
the instant of ionization. 

Concerning the spin asymmetry during ionization, one may 
ask two independent questions: 

1. Does the ionization rate depend on the initial spin state 
of the bound electron? 

2. Will the electron be polarized after tunnel-ionization 
from an unpolarized target? 

Consequently, we can define two spin asymmetry parameters, 
firstly, the tunneling asymmetry parameter 


Let us first consider predictions of s-cc-SFA, which neglects 
the spin dynamics in the bound state. We underline that due 
the latter the results of s-cc-SFA are applicable for an arbitrary 
elliptical polarization 


1. The spin quantization axis is parallel to the laser propagation 
direction 

When the spin quantization axis is along the laser propaga¬ 
tion direction z, the angles determining the quantization axis 
are 0-0 and <p — 0. In this case we obtain for the spin 
flip relative probability corresponding to the maximum of the 
momentum distribution 




+ iy+^- - 
Wt 


(98) 


which measures the asymmetry between the tunneling rates for 
initially polarized states. Here 




dW,. 


(fp 


(99) 


P=Pm 


is the differential ionization rate at the momentum value p 


Pi^ corresponding to the maximal tunneling rate, see Eq.(95 1 , 
and 


Wt = 




-H W H- + w_ 


( 100 ) 


^^(^g/2+p)/2 

^ 1 -(- ^ 2 ^ 2/4 


+ 0(p^). 


(103) 


in the leading order of p s ^J2Ip|c^. On the one hand, the 
spin flip vanishes for a linearly polarized filed for the maximal 
tunneling probability. On the other hand, for a nonvanish¬ 
ing laser’s polarization (^ 0), the spin flip is negligible in 

the weak field regime ^ « 1, ^ 0, whereas s-cc-SFA 

predicts almost complete spin flip for strong fields ^ » 1, 

The tunneling asymmetry parameter and the spin polariza¬ 
tion asymmetry parameter can be calculated up to the nonvan- 
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ishing order of p as 


aiZ (i) _ ^ 3 

A 




1 +^2^2/4 


(104) 

(105) 


The s-cc-SFA calculations predict a rather large spin polariza- 
tion asymmetry because the parameter J{p scales linearly 
with the parameter p, while the tunneling asymmetry parame¬ 
ter is much smaller ~ p^. Nevertheless, the spin asymmetries 
disappear for the linear polarization case. 

We note that on the one hand the relation = p^j indi¬ 
cates that there is no singularity in Eq. ( |104 1 in the weak field 
regime (^ « 1), where p « y « 1. On the o ther hand, the 
relationp^/^ = w/c^(£’fl/£’o) indicates that Eq. (104 1 vanishes 
for low frequencies. When neglecting the terms of the order of 
in the limit ^ 1, one has also to neglect the terms of the 

order of p^, because they are much smaller than the former by 
the factor of y^ 1. 


Spin effects for strong fields f » 1 



k 

-£(77, 

= 0) 

1 

II 

O 


s-SFA 

d-SFA 

s-SFA 

d-SFA 

s-SFA 

d-SFA 



1/2±CpV4 

p2/4 

1/2 

1+ 

1+ 


0 


0 

0 

2p 

a-O p" 

J{p 

0 

0 

0 

0 

2p(l - 20 

p"0 - 20 


TABLE 1. Spin flip and asymmetries in strong field regime ^ » 
1 at the instant of ionization associated to the maximal tunneling 
probability in the leading order of the parameter 1/^. Comparison of 
the s-cc-SFA and d-cc-SFA results for different spin quantization axes 
at the leading order inp = ^/27p/c; f = 0 is for linear and f = 1 is for 
circular polarization of the laser field. 


The spin polarization asymmetry parameter as well as the tun¬ 
neling asymmetry parameter scale as p. The latter is also inde¬ 
pendent from the polarization as well as the intensity parameter 
of the laser field. 


2. The spin quantization axis is perpendicular to the laser 
propagation direction 


In the case when the spin quantization axis is perpendicular 
to the laser propagation direction, one can align it along the 
direction of the electric field or the magnetic field at the in¬ 
stant of ionization rjr, and for the sake of convenience we will 
consider the directions opposite to the fields. 

When the spin quantization axis is along the direction of 
-Eirjr -Q) {9 - 7r/2,0 = 0, see Fig.[^, cf. Eq. (81 1 , the flip 
probability is 


rX-JC (J) _ 
Z ± ~ 



(106) 


which is tiny for nonrelativistic (low charge) ions. Eurthermore, 
the asymmetries are vanishing 


^0. (107) 


The results are independent from both the polarization and in¬ 
tensity parameters at the maximum of the tunneling probability, 
L = 0. 

When the quantization axis is along the direction of -Bitjr - 
0) (6 - njl, <p - njT), cf. Eq. ( [82] l, the spin flip is 


(^¥/4)(1+p) 

1 -H ^2^2/4 


+ f5(p'). 


(108) 


While it is insignificant for weak fields, a complete spin flip 
occures in the strong held regime when 4" 4^ 0, which is similar 
to Eq. ( |103[ ). On the other side, the spin asymmetries can be 
written as 


Jit, = 2 p, 

^■v(b . 1-C¥/4 

'^p -^^ 1 +^ 2 ^ 2/4 


(109) 

( 110 ) 


B. The d-cc-SFA prediction 

In the weak held regime ^ « 1, as the asymmetries as 
well as the spin flip are negligible (smaller by an order of 
magnitude than p^) for any choice of the quantization axes, 
in what follows, we discuss the predictions of the strong held 
limit of d-cc-SEA. Note that in contrast to the s-cc-SEA case, 
the polarization parameter is ^ = {0,1) for the d-cc-SEA case. 


1. The spin quantization axis is parallel to the laser propagation 
direction 


In this case the spin flip in the strong held regime (^ » I) 
in the leading order of p equals to 


<jrZ (d) _ 



( 111 ) 


which is due to the spin fast oscillation in the bound state (see 
Sec. |VIB' for an intuitive description), and it is different than 
the prediction of the s-cc-SEA, whose strong held limit yields 
^z (s) ^ ^ _ jQ^ j puj-thermore, there is a correction 

at the order of p^ for the circular polarization case, see Eig.[^ 
The corresponding asymmetries are 


= (p^ ( 112 ) 

=6)(^-'). (113) 


Eirstly, the polarization asymmetry parameter vanishes for both 
of linear and circular polarizations. The tunneling asymmetry 
parameter also disappears in the linear polarization case, where 
s-cc-SEA and d-cc-SEA agree with each other for strong fields. 
However, the d-cc-SEA result for the tunneling asymmetry 
parameter differs from the s-cc-SEA one in the case of circular 
polarization, as it increases with increasing charge of the ion. 
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2. The spin quantization axis is perpendicular to the laser 
propagation direction 

First of all, when we align the spin quantization axis along 
the -Eirjr - 0)-direction, the spin flip can be calculated as 

(114) 

which is again a consequence of the bound state dynamics, cf. 
Eq. P06| l. However, as in the case of s-cc-SFA, the asymme¬ 
tries vanish 


= = (115) 

In the case where the spin quantization axis is the -B{T]r — 
0)-direction, the spin flip is 

= f (l± yj + ar*)- ( 116 ) 

We first notice that the spin flip vanishes for the linear po¬ 
larization case. For a circularly polarized held the spin flip 
is the same as in the case of the quantization axis along the 
^-direction. 

The asymmetries in the d-cc-SFA case are 


man model spin transitions taking place at an arbitrary instant 
of ionization. 

The complete propagation of the tunnel-ionized state can be 
decomposed into three parts as 


f/(oo, -oo) = t/c(oo, IrWrits + tr, trWgitr, -“) , (119) 


where Ugitr, -oo) describes the propagation of the bound state, 
Urits + tr, tr) the under-the-barrier dynamics, and Uc{°°, tr) the 
continuum dynamics, see F ig.[^ Here, tr is an arbitrary instant 
of ionization, and = i yjllplEQ is the imaginary Keldysh 
time given by Eq. ( |96| l, which defines how long the electron 
has to travel for the width of the tunneling barrier in imaginary 
time and determines the tunneling rate. 

Firstly, the spin resolved continuum evolution can be found 
via the spin resolved prefactor of the Volkov wave function Q 
in the following way. The Volkov bispinors in two different 
times are connected to each other as 


I + -— (l -H ^ ■ O') Ait) ■ a 


2cA 

= licit, t') 


( 120 ) 


I H- (l H- ^ ■ or) Ait') ■ a 

2cA '' ' 


Vj, 


with the identity matrix I. Hence, the continuum time evolution 
operator yields 


^-;'''^-(i-0p'+o(r‘), (117) 

( 118 ) 


In the previous configuration, Eqs. ( |I I2| l and ( I I3| l the both 
asymmetries vanish for a linearly polarized held, however in 
this case, the both have the same nonvanishing value. More¬ 
over, we notice for the circular polarization case that when 
the spin quantization axis is rotated from the z-axis to the - 
Bipr = 0)-direction, the tunneling asymmetry parameter and 
the spin polarization asymmetry parameter interchange their 
roles except a sign difference. Namely, = 0, 

and - - Jit, - p^. 

Comparison of s-cc-SFA and d-cc-SFA for strong fields, 
where the spin effects are not negligible, are presented in Ta¬ 
ble |I] in a compact way. Thus, the standard and dressed SFA 
give different answers for spin asymmetries and the question is 
which result has physical implication. In the next section we 
develop a simpleman model for the spin dynamics in tunnel- 
ionization which allows us to identify the origin of the differ¬ 
ence between s-cc-SFA and d-cc-SFA, and to recognize the 
d-cc-SFA result as the physical relevant one. Note that the 
same conclusion is reached in ll59l from the comparison of 
analytical results with numerical simulations in the case of a 
linearly polarized laser held. 


Ucit, T) - I + ^ (l + * ■ «) [A(f) - Ait')] a, (121) 
which reduces to 

Uci°°, fr) = I - ^ (l + * ■ O') Aitr) a, (122) 

where we use the fact that A(oo) —> 0. Moreover, we will 
consider only the large spin components, therefore, moving 
from the bispinor- to spinor-description 

Uci°°, tr) ~I - [sin(wf^)cr 3 , -H ( COS(cufr)cr;t] , (123) 


3(0 



VI. SIMPLEMAN MODEL EOR THE SPIN DYNAMICS 

All the SFA results calculated in the previous section con¬ 
cerning the spin dynamics can be intuitively inferred in the 
following simpleman model. We will describe by the simple 


FIG. 3. (Color online) The simpleman model for the time evolution 
of the initial bound spin state. The propagation can be split up to 
three parts; the bound state, the under-the-barrier, and the continuum 
dynamics. The standard SFA (red dashed curve) neglects the former, 
whereas the dressed SFA (blue solid curve) takes the laser field into 
account in the bound state dynamics. 
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where we further set /I = 1 for the spinor-description consis¬ 
tency. 

In order to model the time evolution for the under-the-barrier 
dynamics, we consider a triangular barrier depicted in Fig. 
and apply WKB approximation. The WKB propagator can be 
mimicked by the following spin resolved propagator, which 
is derived via including the spin interaction Hamiltonian into 
Eq. 


Urits) - exp <1 dt 

with e± being the energy of the tunneling spin-1/2 particle, 
which is different in the standard and dressed SEA. 

Since the bound state dynamics varies, we will investigate 
it in a more detail way for the standard and dressed SEA, 
respectively. 


q(ty 


■ S± -t- 


o- ■ B{t) 
2c 


, (124) 


A. Standard SFA 

In the standard SEA, the influence of the laser field on the 
bound state is neglected. Due to the latter, all the results which 
we derive are valid also for an arbitrary elliptical polarization. 

The evolution of the bound state between two different times 
in this case is given by the usual propagator as 

Usit, t') - exp [-/e±(f - f')] = exp - f')] , (125) 

with the energy of the tunneling bound state particle s± = -Ip, 
which leads to a trivial phase. Accordingly, as there is no 
Zeeman splitting of the energy in the bound state, the under- 
the-barrier propagator reads 


Urits) = exp 


Jtr 2C 


(126) 


where the tunneling probability amplitude yfWiTr) arises from 
the two hrst terms in Eq. \\2A\ , which depends on the arbitrary 
instant of ionization (or the electron’s final momentum). Eor 
example. 


yW(0) - exp 



(127) 


for the instant of ionization that corresponds to the maximal 
tunneling probability. Eurthermore, the spin resolved term can 
be written up to the order O(t^) as 


rdt^ 

Jc 2c 


* (^<^x sin(wfr) - (Ty cos(wfr)) , 


2c 


(128) 

(^CTy COS(OJtr) - sin(wfr)) . (129) 


tion axis s can be written as 


+ (oo); s) = ^JW(tr) * ^ (sin(mf^)crj, -H ( cos(ojtr)o-x) 


X exp 


(cTj, cos(wfr) - ^CTx sin(mfr)) 


l±;5), 

(130) 


and the transition amplitude can be given by 

= (s'l [l - (sin(«fr)cry -i- ^ cos{Lotr)o-x) 


(131) 


X exp - (^cTy cos(wfr) - ^CTx sin(mfr)) 




where we omit the tunneling probability amplitude for the sake 
of simplicity. 

Let us hrst consider the case when the spin quantization axis 
is chosen along the z-direction. The transition amplitudes up 
to the leading order in p can be found as 

= 1 - 1 - \±2^ + - 1) sin(2mf;.)j , (132a) 

i [(+^ - ^p) sin(«fr) - + p) cos(mfr)] . 

(132b) 


Eor the spin Hip, we derive 




_4(2 + ^^p)_ 

8 -H (1 H- - (1 - /' 2 )f 2 cos(2wf,) 


. (133) 


Eor the instant associated to the maximal tunneling probability, 
fr = 0 we have 




^^(^^ + 2p) 

4 -H ^2^2 


(134) 


which coincides with the s-cc-SEA result given by Eq. (103 i. 
If we go further in the next orders of p, the spin Hip at = 0 
can be written as 


jrz(s) 


(0)^ 


4 + ^2^2 • 


(135) 


Eq. ( |135^ shows that both the spin precession in continuum (the 
term ~ ^) and the spin dependent tunneling probability (the 
term ~ p) are contributed for the spin flip effect. In the case 
of linear polarization, ^ — 0, Eq. (135 i gives Tz/ *^*^(0) p^/4, 

indicating spin flip only due to the tunneling. However, this is 
specific only for = 0. Eor this reason one can introduce the 
mean value of the spin flip over the laser’s period To, which 
yields 






^(^ + /-2^ + 4^p) 

8 + (l+/'2)f2 


(136) 


ZTT 


(jjln 


dt VT,. 


As a consequence, the final spin state in terms of the initial 
spin state |±; s) = | ± (-oo); s) with an arbitrary spin quantiza- 


with 


(137) 
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The average spin flip in some cases can be different from the 
instantaneous value. For instance, in ^ = 0 case, Eq. ( |136| l 
provides the strong held asymptotic, f » 1, = 1 in 

contrast to p^/4. 

At the same order the spin asymmetries can be given by 

(138) 

Wtf -) =__ cnQ-) 

^ 8 + (l+^2)^2_(l_^2)^2cos(2wf,)' ^ ^ 

The tunneling asymmetry parameter is negligible for all possi¬ 
ble intensities as well as laser’s polarization similar to Eq. ( |104| i. 
Further, it is independent from the instant of ionization. The po¬ 
larization asymmetry parameter disappears for both of the weak 
held and the strong held regimes. Nonetheless, the polarization 
asymmetry parameter depends on the ionization moment for 
intermediate fields. For the maximal tunneling probability, it 
reads 


T7#74 • 

which agrees with Eq. ( |105| l, whereas its mean value over the 
laser’s period can be calculated as 




1+(1+^2)^2/8’ 


(141) 


7 ^5^ 

which has the same qualitative behavior as Jip (0). 

In a similar way, we can deduce by the simpleman model 
the spin flip as well as the spin asymmetries for different spin 
quantization axes. Eor example, when we align the quanti¬ 
zation axis along the y-direction, we provide the following 
expressions up to the order 0(p) 


y (j) ^ 2 ^ 2^2 cos(ajt,)(cos(ojt,) + p) 

^ 8 -H (1 H- ^ 2)^2 _ (1 _ ^ 2)^2 cos(2a>tr) ’ 

^*Vr) = 2pCOS(a>tr) , 

= 2pCOS(a>tr) 


(142) 

(143) 

(144) 


X 


4^2^2 

8 -H (1 H- ^2)^2 _ (1 _ ^2)^2 cos( 2 co>tr) 


Eor the linear polarization case, the spin flip vanishes and 
the both asymmetries simplify to - 

2p cos(wfr). The spin asymmetries as well as the spin flip for a 
nonvanishing polarization depend on the instant of ionization 
even for strong fields. While we derive 


Fl ^'^(0) 


(l±p)^¥ 

4 + ^ 2^2 ’ 

2p, 

2p(4-^2^2) 
4 ^2^2 


(145) 

(146) 

(147) 


Eqs. (1081-( 110 1 , their mean value can be calculated as 


riel 

irpy (b\ _ S _ 

^ ’ 8 + (l+^2)^2’ 

= 0 , 


(148) 

(149) 

(150) 


One should note that the time averaging destroy the asymme¬ 
tries. The spin flip qualitatively remains the same at the time 
averaging, but decreases two times. 

If the quantization axis is along the .r-direction, at the leading 
order of p, the spin flip and the spin asymmetries can be written 
as 


^ sin(fa;t,)(sin(cut,) + ^p) 

* ^ 4 + ^2 (^2 cos2(wf,) + COs2(wf,)) ’ ^ ^ 

Ji] = -2^p sin(wf,), (152) 

^ _ 2 ^p sin(wf,) (153) 

„ r 4^2 

8 -H (1 H- ^ 2)^2 _ (1 _ ^ 2)^2 cos(2wf,) ' 

All the spin effects depend on the instant of ionization, and all 
of them disappear for the maximum tunneling probability at 
the leading order in p. If we further go the next order, the spin 
flip for the maximal tunneling probability can be given by 




(154) 


which is consistent with Eq. (1061. The mean values over the 


laser’s period, on the other hand, can be given by 


^ ^ 8 + (l+^2)^2’ 

(155) 


(156) 


(157) 


In this case, the average value of the spin flip is quite different 
than the instantaneous value. Thus {F^ 1 at ^ » 1. 

We notice that our simpleman model captures the main fea¬ 
tures of the spin dynamics and is able to reproduce the results 
of the SEA. The averaging of the spin flip and asymmetries 
over the ionization time in some cases modifies significantly 
the result. 


B. Dressed SFA 

In the dressed SEA, the scenario is changing because we 
take into account the bound state evolution in the laser field. 
The bound state propagation is not a trivial phase in this case, 
but includes the spin precession in the bonnd state. Accord¬ 
ingly, the Hamiltonian for the bound state propagation can be 
represented as 


H ^ Hoi + 


a- ■ B{t) 


for the maximal tunneling probability, and we notice that 
the simpleman resnlts captnre the corresponding SFA resnlts. 


2c 


(158) 
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FIG. 4. Intuitive description of the spin dynamics under-the-barrier 
can be modeled with a triangular potential barrier. In the standard SFA 
case, the spin states of the tunneling particle enter the barrier with the 
same energy —Ip and are split during the under-the-barrier dynamics, 
whereas in the dressed SFA, first the energy of the spin state is split 
due to the Zeeman effect in the bound state, but later it is compensated 
by the energy splitting in the under-the-barrier dynamics. 


First of all, the transition amplitude is a function of a-y, as 
a consequence there cannot exist a spin flip when the spin 
quantization axis is the y-direction, which explains the derived 
result ( |116| l. Furthermore, due to the symmetry reason, the 
spin flip for the case when the spin quantization axis is the 
z-direction is the same as when it is the .r-direction, and they 
can be calculated as 

* ^ 8-t-^^(1 - cos(2wff)) 

(163) 

At tr = 0, ^‘'^0) = 0. If tp + 0, Ti"" ~ for weak 

fields, and it is negligible. Flowever, its strong held limit can 
be given by 

= cos^ sin(wf,)) . (164) 


where Hq is the usual atomic Hamiltonian of the electron in 
the Coulomb held of the core. We would like underline that 
once the solution of the Schrodinger equation for the Hamilto¬ 
nian ( |158[ ) is known, the developed simpleman model is valid 
for any polarization, and accordingly we discuss the linear and 
circular polarization cases separately for the dressed SFA. 

In addition to the trivial phase ( |I25 I, in this scenario the 
bound state propagator for a linearly polarized held can be 
calculated as 




Usitr,-^) - exp -/ 


= cos sin(wfr)j I -I- icTy sin sin(wC)j . 

(160) 


O-yByif) 

2c 


(159) 


Furthermore, due to the Zeeman splitting, the energy of the 
tunneling bound spin states whose quantization axis is along 
y-direction will read s± = -Ip + Eq cos(cufr)/(2c) for a certain 
instant of ionization C, where s+(_) is for spin-up (down) state. 
Consequently, as the spin states having different quantization 
axes can always be expanded by the basis vectors |+;y), the 
Zeeman splitting and the spin interaction under-the-barrier 
cancel each other, and the propagator is simply given by the 
tunneling probability amplitude 

uAts)^ ( 161 ) 


We conclude that the tunneling probability is spin independent 
in the dressed SFA at the leading order in p. Therefore, the 
spin asymmetries disappear at this order [the contribution of 
higher orders are discussed in the next subsection]. 

After including the continuum propagation for ^ = 0, the 
transition amplitude in the case of the dressed SFA for linearly 
polarized held becomes 


= x<i 


'l(l- sin(mC)o-3,) 


(162) 


^cos sin(wf^)j I + icTy ^2 ®l*^(^C)jj k) ■ 


As we discussed in Sec. Ill the bound spin highly oscil¬ 
lates between up and down states in a very short time interval 
for strong fields. As a result, it is physically more correct to 
present the mean value of the spin flip instead of its instanta¬ 
neous value at C = 0. Furthermore, one can obtain the mean 
value corresponding to the maximal tunneling probability by 
averaging over a period T fulfilling the condition (|74|), and it is 


(d) 


)t - 


1 

2 ’ 


(165) 


which agrees with Eq. ([TTT), and Eq. ( [TT^ as well as with 
Ref. 1591 . Note that the mean value over the laser period Tq 
needs a numerical calculation. 

Eor the circular polarization case, we first notice that solv¬ 
ing the coupled differential equations ( |54| l corresponds to 
solving the Schrodinger equation for the Hamiltonian H — 
<T ■ B{t)l{2cl6). Eurthermore, the corresponding propagator 
can be written as 


'C++(f,) C-+(f,)' 
C-+(f,) C-(f,), 


(166) 


with the coefficients ( |66| l. As a consequence, the bound state 
propagator for a circularly polarized held can be obtained with 
the replacement of ^ > 1 in the coefficients ( 661 , which reads 

1 




^^2 + (i+ 


(167) 


(1 + 


iojti 
i+y/i+£- 


1 -V 1 +? 


i^e 


1+Vi+F 


(1 -H 


i-VFF 


While the under-the-barrier propagator is the same as Eq. ( |161| l, 
the continuum propagator for the circular polarization case 
reads 


Uc(°°, tp) - I - |^sin(wfr)(Ty + COS(wfr)cr;cj . (168) 
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Thereby, we can give an intuitive description of the spin effects 
for a circularly polarized field. For instance, when the spin 
quantization axis is along the z-direction, the spin flip is 




^2 [4 - 4(1 + E)Cos (Wh(l - ^)) + (1 + Ef\ 


2(4+ ^2) ( 52 +5) 


(169) 


with E = For a weak field the spin flip is negligible, 

whereas in the strong held, f » 1 we derive 


n 


z(d) 


( 0 ) 4 ’ 


(170) 


for the instantaneous value at tr - 0. The mean value of the 

1 


spin flip over the period T (74 1 , 




(171) 


which coincides with its instantaneous value. This result is a 
consequence of this particular choice of quantization axis as 
we discussed in Sec. Ill A In a similar way, the strong field 
limit of the spin flip for the other spin quantization axes can be 
derived as 


1 - 


- - 


While at h = 0, we provide 


cos [«fr(l - S')] - COS [«h(3 - S)] 



T± ^“'\0) = 0, 


their mean values over the period T are 




1 


(172) 


(173) 


(174) 

(175) 

(176) 

(177) 


which are consistent with Eq. ( |1 16| l, and Eq. (114 1 . 

Thus, the simpleman model clearly shows that the spin dy¬ 
namics in the bound state is very important for developing 
of spin effects during tunnel-ionization. This dynamics is in¬ 
cluded only in the dressed SEA and, therefore, the physically 
correct results are those predicted by the dressed SEA. 


C. Origin of the spin asymmetries in the dressed SFA 


If we compare the results developed in the previous sec¬ 
tions VIA and VIB we notice that the origin of the spin 


asymmetries in the standard SEA is the spin interaction under- 
the-barrier, which leads to different tunneling probabilities for 


spin-up and spin-down states. Since the Zeeman splitting in 
the bound state compensates this spin interaction in the dressed 
SEA, there is no spin asymmetries at the leading order of p. 
However, the magnetic field in the rest frame of the electron 
is slightly different in the bound state and under-the-barrier, 
which have an impact on the Zeeman splitting compensation 
and can lead to a spin asymmetry. This can be described by 
improving our simpleman model in the following way. 

If we go further in the Eoldy-Wouthuysen expansion of the 
Hamiltonian, we can write down the propagator for the under- 
the-barrier dynamics as 


Urits) 

- exp 


(. rtr- 


(178) 


dt 


q(tf 


„ , o-.B(f) , o-{Emq{t)) 


}■ 


where the last term in the exponent of Eq. ( |178| l describes the 
effect of the magnetic field in the rest frame of the electron 
during tunneling. Note that as the momentum vanishes for 
the bound state, i.e., {p)b - 0, we do not need to modify the 
energy of the tunneling bound spin states. Moreover, the rest 
frame effect depends on the final momentum distribution via 
the term E{t) x q{t). Thereby, in order to elaborate the spin 
asymmetries in the dressed SEA, we consider the instant of 
ionization associated to the maximal tunneling probability. The 
rest frame effect term, then, can be calculated up to the order 
Oii]) as 


‘i - 


where we have further used q^(t) - -2Ipl{3c) via Eq. ( [24| . As 
a result, the under-the-barrier propagator yields 

Urds) = 


.o-yEoq^(0)t, 

4c2 

O-yP^ 

12 ’ 


(179) 


Ub { Q , -«)) = ( 1-^)1 + 


M, 


= <i'l 




(1 - ^)i ■ 


^e+{\ + yji+e)^ 


(181) 


for the maximal tunneling probability. 

Combining the bound state propagator for linear ( |159| l and 
circular polarizations ( |167| l together, we derive 

^[(1+ VrTF)l + f^cr,] 

with ^ = {0,1). Including the continuum propagator 
t/c(oo,0) = I - lyo-x, 
the transition amplitude can be written as 


(182) 


(183) 
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where an expansion over a small parameter p is used. Then, 
in the linear polarization case, ^ - 0 the transition amplitude 
yields 

Since the amplitude depends only on o-y, there cannot exist the 
spin asymmetries when the spin quantization axis is along the 
z or ;ic-directions. This is the underlying reason of Eqs. ( |1 12 1 
and ( |113| l as well as Eq. ( |115| l. However, when the spin quanti¬ 
zation direction is chosen in the y-direction, we derive 


ts i yjUplEoil - 5p^/72) as well as the effect of the bi¬ 
spinor description can improve the accuracy of the results of 
the simpleman model. 

Next, we have neglected the spin-orbit interaction during 
the electron motion in the continuum under the action of the 
laser and Coulomb fields. In the paper II, we have identified 
the spin-orbit coupling term in the continuum wave function. 
Including this term in the continuum propagator will result in 
a new term for the matrix element describing spin transitions 

1 r“ 

— (188) 


tziy (^)_ fip ^ 


(185) 


which agrees with Eqs. (117 1 and ( 1 18| l up to 1/3 prefactor. 

In the case of a circularly polarized held, the corresponding 
spin asymmetries can be written as 




rZ (^) _ 


3^/TTJ■ 




y (d) 


3^/TTJ■ 


:P^ , ^ 


3 (Tiy {d) _ 

:P , - 


4^ 


3iA+ey 

4_^2 
3(4+ ^2/ 


= 0 , 




X (d) _ 


= 0, 


(186a) 

(186b) 

(186c) 


and their strong held limit, ^ » 1 is 

„3 

^z(d)^P_^ (187a) 

yi; = 0, , (187b) 

=0. (187c) 


Similar to the linear polarization case they agree with the d-cc- 
SFA results up to the 1 /3 factor (see Table|I]i. 

As a conclusion, the simpleman model describes correctly 
the qualitative behavior of the asymmetry parameters. Here 
we would like stress that we divide the total propagator into 
three parts: the bound, under-the-barrier, and continuum parts. 
While the bound and continuum propagators are unitary opera¬ 
tors, the tunneling asymmetry parameter JTI, vanishes for these 
propagations. The latter follows from the fact that since the 
conservation of the probability implies W+^+ + = 1 as 

well as W _>+ -I- VT—>_ = 1, with = Ks'|t/(f,t')|s)l^, by 

dehnition essentially disappears. Nonetheless, the under- 
the-barrier evolution operator is not a unitary operator as it 
determines the tunneling rate. Therefore, the asymmetries are 
originated from the tunneling step. 


D. Further contributions to spin effects 

Finally, we would like to discuss further contributions to 
the spin effects that have been neglected in our treatment. 
First of all, we would like to emphasize that the following 
relativistic corrections, T —» 1 - 16, <5^1- p^ 13, and 


This will modify not only the simpleman model but also the 
SEA calculations. 

Finally, there could be a relation of spin effects to the tun¬ 
neling delay time ll8^ . which needs investigation beyond the 
quasiclassical description of tunneling. For instance, in the 
simpleman model the time evolution of the wave function will 
be modihed as 


U{oa,-oo) - Uc{°°,tr+T)UT(tr + ts+T,tr)UB{tr,-°°) , (189) 

if we introduce a real and positive time t associated to the 
time spent under-the-barrier, which would lead to an additional 
precession and spin effect modihcations. 

The mentioned modifications for the spin transition proba¬ 
bilities will be discussed elsewhere. 


VII. EXPERIMENTAL OBSERVABILITY 
A. Detection of photoelectrons 

The measurement of the photoelectron spin flip during 
tunnel-ionization requires an initially polarized atomic target 
and a detection of the photoelectron spin polarization. The 
photoelectron spin polarization can be measured using Mott 
polarimetry ll83H85]l . The latter is based on the left-right asym¬ 
metry of Mott scattering cross-section on a high-A- target. It 
depends on the electron spin polarization !Px, transverse to the 
scattering plane; dcr — derail + !Pj_5(0)], where dcr and do-Q 
are the spin resolved and spin-averaged cross-sections of Mott 
scattering, respectively, S (0) is the so-called Sherman function, 
the maximum of which is approximately 0.5 in the case of a 
gold target at electron energies of the order of megaelectron- 
volts HH, and ~ 10 2 in the case of zinc or lead targets at 
electron energies of order of hundred megaelectronvolts Il84l . 
For instance, when the spin of the ionized electron is polar¬ 
ized in the laser propagation direction in a circularly polarized 
laser field, the electron spin polarization in the final state is 
1/2 for ^ » 1, according to Table |l] and see Fig. This 
means that the spin in the final state will be oriented trans¬ 
versely with respect to the propagation direction. Therefore, 
~ 1 is possible, because in the relativistic regime the elec¬ 
tron final momentum is mostly along the laser propagation 
direction. At f ~ 10 (the laser intensity of the order of 102** 
W/cm2), 'Pj^S (0)max ~ 10 ^ 2 ^ and the relative error of the signal 

(~ 1/ ■JNs \ with the total scattering events A^') should be 
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FIG. 5. (Color online) The spin flip associated to the maximal tunnel¬ 
ing probability for a circularly polarized field when the spin quanti¬ 
zation axis is the laser’s propagation direction within d-cc-SFA. The 
blue solid and the red dashed lines represent and re¬ 

spectively. The applied parameters are tu = 0.05, and EqIE^, = 1/30. 


smaller than 10 ^ to distinguish the electron polarization, i.e., 
= VshotNs^^ > 10"^, where Vshot is the number of laser shots, 
and is the number of scattering events per the laser shot. 
The latter is determined by the ionization (Wi) and the Mott 
scattering (Wm) probabilities ~ Wi x Wm x (dOln), where 
the effective interval of the scattering angle (where S (0) ~ 10“^ 
f85ll ) is 66 ~ 5°. Estimating the Mott scattering probability via 
Wm ~ ~ '■h® Bohr radius gb, the target 

length ^ = 0.1 /tm, the Lorentz factor yi ^ - 100, k - 30, 

and the density pQ - 2 ■ 10^^ cm“^, each electron will produce 
a scattering event but the hindering multiple scattering will 
be avoided. The ionization probability per laser shot can be 
estimated as Wi ~ (at a fixed EjEa - 1/10, the laser 

pulse duration of 100 fs) ll66l . Then the scattering events per 
laser shot is Wm x {56jn) x VT,- ~ 3 x 10“^. The number 

of laser shots required for the necessary signal resolution is 
Vshot > nP/N s ~ 3 X 10^, which can be realized even with 1 
Hz laser system. 

For typical experimental parameters, e.g., ionization of hy¬ 
drogenlike Ne®^-ions in a strong infrared laser field with an 
intensity of 10^*^ W/cm^, the d-cc-SFA predicts a spin-flip 
relative probability of about 0.1. 


B. Detection of ions 


being the orbital angular momentum and the spin of the ion and 
electron, respectively, andziy, ,, = J/g- Ji,e, i^e - S ^ - S i^e 
and n - AJe + A Ji. Taking into account that the change of the 
electron angular momentum during above-threshold ionization, 
which is 


AL, 


In + 2Un 


U) 


(192) 


see 112, the change of the electron spin can be measured via 
the ion angular momentum 


AS = -AJi. 


(193) 


Let us employ a circularly polarized laser held and polarize the 
initial electron spin along the laser propagation direction (the 
ion nuclear spin is assumed to be vanishing Si - 0). Even when 
the electron spin is not changed during the ionization AS e = 0, 
the ion carries out a small part of the angular momentum 
provided by n absorbed photons. In particular, with a right 
circular polarization, the ion gets an orbital angular momentum 


Lf = Lo = 


in Ip 2 Up 

M a> 


(194) 


where m and M are the electron and ion mass, respectively. 
The latter follows from the fact that the ion distance with re¬ 
spect to the center-of-mass of the ion-electron system m/M 
times smaller than the electron distance. With a left circular 
polarization, the ion acquires an angular momentum = -Lq- 
When the spin flip happens, - Lq - \ and - -Lq - 1, 
respectively. Therefore, the spin flip will be indicated by a non¬ 
vanishing signal for the difference in the ion angular distribu¬ 
tion with the left and right polarized laser held. The mentioned 
signal for the difference in the ion angular distribution can be 
estimated as 


‘S~||Tl„+i,Lo+iI"-|1'a,-i,l„-iI"I (195) 

~ I sin^'^+'l 6 - sin^'^''^*’ 6\ ~ 26, 

where T; are the spherical harmonics, 6 is the ion scattering 
angle with respect to the laser propagation direction, and 6 = 
cot^ 6* 1. Then, the required number of laser hots is Vshot ~ 

10 ® at A ~ 10“^ (the ions are observed in the transverse to the 
laser propagation direction within an angle of 20 mrad). 


The electron spin flip can also be revealed via a measurement 
of ion parameters. The angular momentum change of the ion 
during ionization AJi can be related to the electron spin change 
AS , using the angular momentum and the energy conservation 
laws in the case of circular polarization 

n + Ji + Je^J' + J'„ (190) 

na)^Ip+2Up, (191) 

where n is the number of the absorbed photons at ionization. 
Up - 12 is the ponderomotive potential, Ji e and 7'^ are the 

ion and the electron total angular momentum before and after 
the ionization, respectively, Ji e - Li e + S with and S i^e 


VIII. CONCLUSION 

In this paper we have investigated the spin effects during 
strong held ionization in the tunneling regime for linearly as 
well as circularly polarized laser fields. The spin effects are 
fully described by two types of spin asymmetries, tunneling 
asymmetry and polarization asymmetry, as well as by the spin 
flip probability. The tunneling asymmetry describes the asym¬ 
metry in ionization of the initially polarized target, while the 
polarization asymmetry describes the degree of the photoelec¬ 
tron polarization from an unpolarized target. The spin resolved 
differential ionization rates as well as the spin asymmetries 
are derived for the maximal tunneling probability. For this 
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purpose, the final momentum distribution of the tunnel-ionized 
electron is calculated. The results are further generalized to an 
arbitrary spin quantization axis. 

Two versions of the Coulomb corrected SFA have been con¬ 
sidered. While in s-cc-SFA the influence of the laser field on 
the bound state is completely neglected, the latter is properly 
treated in d-cc-SFA. Therefore, we have discussed the spin 
dynamics in the bound state driven by linearly as well as cir¬ 
cularly polarized fields. The physically relevant predictions 
for the spin effects is given by the dressed SFA. Generally, the 
spin effects calculated with s-cc-SFA are overestimated. 

A simpleman model for spin effects is developed which 
transparently shows how the spin effects arise in three steps; 
spin precession in the bound state, spin rotation during tunnel¬ 
ing and further spin precession in the continuum. 

The parameter ^ determines the field regime, and in the weak 
field regime ^ 1 the spin flip as well as the spin asymmetries 

are negligible. The spin effects are considerable for strong 
fields ^ » 1 , where the spin flip probability is about 1 / 2 . 
The spin flip probability is mainly determined by the bound- 
state dynamics, however the latter has no effect on the spin 
asymmetries. Indeed, the spin asymmetries, which increases 
with increasing the charge of atomic core, are a consequence of 
the effect of the magnetic field in the rest frame of the electron 
during tunneling as indicated by the developed simpleman 
model. 

The most favorable for experimental observation of spin 
effects is the spin flip by using moderate highly charged ions 
with a charge of the order of a: ~ 20 and a laser field with an 
intensity of / ~ 10^^ W/cm^. 
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IX. APPENDIX 


A. Classical equations of motion of a charged particle in a 
plane wave 


e\k - e 2 k - €162 - 0. Then the field strength tensor reads 

+ (197) 

where e^’' = e^'' - k'^e^ - k''^, and dot denotes the 

derivative with respect to the phase 77 . 

The classical equations of motion are obtained by solving 
the Lorentz force law, which can be written in the proper time 
parametrization as 

^ = (198) 

UT C 

The relation k^^^ - kfjel^'' - 0 implies that kp is a constant of 
motion. Using further the following identity 


dpi^ ^ dpf‘ drj ^ 
dr dp dr 


(199) 


with A-kpjoj, the Lorentz force law in the phase parametriza¬ 
tion yields 


//=— FAp. (200) 

c Aojc 

Moreover, the last term in the above equation can also be given 
by 




dp 


dp 


( 201 ) 


where we used the contraction of Eq. (200 1 with A^ after the 
second equality sign. Then, Eq. (200 1 can be written as 

1 




( 202 ) 


dp\c Acjc 
whose solution in terms of the initial phase pi can be found as 
Ai^(p) - Ai^(pi) 


Ain) = Aim) + 


k^‘ 




A\p)-A^rjd 

2c 


iAyip) - Ayipi)) , 

(203) 


where we used the following relation 

Af,ip)Ain) = A^{rf)Aini) + -A^ip) iA'^ip) - A>^ipi)) . (204) 

c 

In the velocity gauge A^ = (0, cA), the three-momentum and 
the energy, then, yield 


The most general gauge potential for a plane wave can be 
given by 


A'' = e^/i(?7) + 


(196) 


pin) = pirn) + Aip) - A(pi) 

^ j . (^Aip) - Aipi)) , (205) 

sin) = sipi) + J ^pipi) + j . (^Aip) - Aipi)). 

(206) 


with the phase p - kxlca and the wave vector k^ iA - 0). Here 
4 , 4 polarization vectors such that - -1 and 


2 
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